Cooperative crosslink (un) binding in slowly driven bundles of semiflexible filaments 
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Combining simulations and theory I study the interplay between bundle elastic degrees of freedom 
and crosslink binding propensity. By slowly driving bundles into a deformed configuration, and 
depending on the mechanical stiffness of the crosslinking agent, the binding affinity is shown to 
display a sudden and discontinuous drop. This indicates a cooperative unbinding process that 
involves the crossing of a free-energy barrier. Choosing the proper crosslinker, therefore, not only 
allows to change the composite elastic properties of the bundle, but also the relevant time-scales, 
which can be tuned from the single crosslink binding rate to the much longer escape time over the 
free-energy barrier. 
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Filamentous biopolymcrs, like F-actin or microtubules, 
play key roles in many cellular and extra-cellular pro- 
cesses. A variety of different crosslinking proteins are 
used in order to assemble these filaments into higher or- 
der structures, like bundles or networks. Permanently 
crosslinkcd reconstituted F-actin networks have been 
used as simple model systems to advance the understand- 
ing of the mechanics of the more complex physiological 
systems Similarly, experiments with actin bundles 

allowed to assess the influence of the crosslink stiffness 
on the bundle mechanical properties jij. A theoretical 
description was developed that characterizes bundle me- 
chanics starting from the assumption that crosslinks are 
permanently bound Recent experiments 0, 0] in- 
dicate, however, that low frequency rhcology of actin 
networks can only be understood by properly account- 
ing for crosslink binding and unbinding processes Q. 
This points the way towards investigating the possible 
coupling between filament elasticity and crosslink bind- 
ing. Several theoretical studies dealt with the equilibrium 
phase behavior of crosslinked networks [§] . Much less is 



known for the case of filament bundles. Refs. (lfj, [ll| 
described an equilibrium unbundling transition in terms 
of a competition between the entropy of filament bend- 
ing fluctuations and the attractive force mediated by 
crosslink binding. Similarly, Ref. 12| invoked filament 
twist elasticity to explain the apparent thermodynamic 
stability of bundles with a characteristic radius 13, 14{ . 



Going beyond unconstrained thermal equilibrium, 
I study here the nonlinear response of a revcrsibly 
crosslinkcd filament bundle to a driving force or defor- 
mation. In the crowded environment of a cell F-actin 
bundles are likely to be constantly exposed to external 
forces and obstacles that bend, compress or otherwise 
deform it. Little is known about how the bundle, and 
in particular the crosslinks inside the bundle, react to 
such a perturbation. I will show that bundle deforma- 
tion leads to crosslink unbinding processes that crucially 
depend on the stiffness of the crosslinking protein. In 



particular, bundles with stiff crosslinks display a cooper- 
ative, discontinuous transition from a highly to a weakly 
crosslinked state. 

Consider a bundle with Nf filaments in two spatial 
dimensions. The bundle is slowly driven into a bent 
state characterized by the tangent angle 8(s) at arclcngth 
position s = . . . L along the bundle backbone (see 
Fig. [1]) . The driving is assumed slow enough such that 
the crosslink binding degrees of freedom can equilibrate 
under the constraint of the given bending amplitude. 
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FIG. 1: Illustration of the bundle geometry. The bundle is 
bent by an arclength-dependent tangent angle 8(s). Filaments 
are discretized into beads, which can move along the bundle 
axis (as quantified by the displacement m a ) but not in lateral 
direction. In addition, each bead can host a crosslink that 
connects it to its neighboring site on the adjacent filament 
(vertical solid lines: bound crosslink; dashed lines: unbound 
crosslink). 



As a consequence of bundle bending, filaments have to 
slip relative to each other, thus bringing the crosslinking 
sites out of registry and leading to crosslink deformation. 
Denote by Ui a the axial slip motion of filament i at ar- 
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FIG. 2: Average crosslink occupation as function of the control parameters fi, k x , c and Nf using N x = 100 and pk s b 2 = 100. 
(left) Undeformed bundle (9 = 0) for different crosslink stiffness k x /k s = 10 _2 ,...,10 3 (from right to left). The solid lines 
are the solutions of the one-crosslink model, Eq. |(3J. (center, right) As function of bundle deformation for different crosslink 
stiffness, k x /k s = 10 _J . . . 10 -2 (center) and for different filament numbers Nf = 2... 8 (right). For the largest crosslink 
stiffness (black circles, center) the simulations display considerable hysteresis, indicating a strongly first-order transition, where 
a free-energy barrier does not permit to sample efficiently. The lines are obtained without fit parameters from the theoretical 
model, Eq. Q. 



clength position s a (Fig. [TJ, then the crosslink (shear) 
deformation can be taken as A^ Q = itj+io, — Ui a + bd(s a ), 
where b is the lateral distance between the two filaments 
that are connected by the crosslink. In the continuum 
limit this reduces to the well known expression for the 
shear strain, du x /dy + du y /dx, with the axial (u x ) and 
the lateral (u y ~ J 6) components of the displacement 
field. Bringing the two crosslinking sites back into reg- 
istry (A — > 0) is therefore only possible if one of the 
filaments stretches out farther than its connected part- 
ner, Ui + \ a = Ui a + b8(s a ), in order to compensate for the 
bending induced mismatch, b0. This leads to the com- 
petition between filament stretching and crosslink shear- 
ing as a fundamental mechanism governing bundle me- 
chanics Filament stretching can only be relieved at 
the cost of crosslink shearing; alternatively, if crosslinks 
remain unsheared, the filaments automatically have to 
stretch out. 

By discrctizing the filaments at the sites of the 
crosslinks one can write the stretching and the shearing 
energy as a sum over all filaments and all crosslinking 
sites as 
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Filaments are thus characterized by a stretching stiffness, 
k s that constrains the axial motion of two successive sites 
along a given filament i. Crosslinks with stiffness k x cou- 
ple two filaments laterally, but at the same arclength po- 
sition, a. This energy expression, which has been given 
previously [!, [Hj], simply corresponds to a discretized 
version of anisotropic continuum elasticity fljj ]. The 
important new ingredient is the possibility of crosslink 
(un)binding, which is accounted for by introducing the 



occupation variables, m a = 0, 1, corresponding to the 
unbound and the bound state, respectively. 

Assuming 9(s) to be constant on the time-scales of in- 
terest, the bending energy of the bundle needs not be 
considered explicitly, as it contributes an unimportant 
constant. The degrees of freedom of the bundle are then 
the axial displacements Ui a of the sites, as well as the oc- 
cupation with crosslinks, I neglect the possibility of 
individual filaments performing lateral motion transverse 
to the bundle axis. It has been shown 13, 11 1 that at suf- 
ficiently low crosslink concentration the entropy gained 
by these fluctuations can drive an unbundling transition. 
Being primarily interested in highly crosslinked bundles 
away from the unbundling transition, these degrees of 
freedom can be assumed to be frozen out 
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This ends the model definition. We will now proceed to 
calculate the average crosslink occupation (n) as a func- 
tion of the control parameters. Combining Metropolis 
Montc-Carlo simulations with theoretical analysis it will 
be possible to elucidate the complex interrelationship be- 
tween crosslink binding and bundle elasticity. 

For the simulations, units are chosen such that b = 
k s = 1. The binding variables nt a enter Eq. linearly, 



so they can be traced out exactly [ll|. As a result one 
obtains an effective crosslink potential that only depends 
on the axial displacements Ui a 
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The depth of the potential is given by j3V(0) = log(l + 
e - ^) and corresponds to the standard free energy of a 
single crosslink with chemical potential \x. The potential 
range is set by the crosslink stiffness, A max ~ 1/ J j3k x . 

Fig. [2k displays the resulting crosslink occupation (n) 
for an undeformed bundle (9 = 0) as a function of fi 
and k x ■ It comes as no surprise that crosslinks tend to 
unbind by either reducing the potential depth (/i /*) or 
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decreasing its range (k x /•). Accordingly, already a sim- 
ple toy-bundle with only one crosslink and two it-degrees 
of freedom allows to reproduce the simulated curves fairly 
well. The Hamiltonian of this simple model reads 

H = y («? + u 2 2 ) + ^ ( Ul - u 2 + b6f+ fin , (4) 

and one obtains for the crosslink occupation (8 = 0) 

-l 



(n) = l + e^Vl + 



2k, 



(5) 



which is shown together with the Monte-Carlo data in 
Fig. 12k. The decreasing crosslink occupation is driven 
by the interplay of binding energy, as characterized by 
/i, and the configurational entropy stored in the axial 
displacement modes. Notably, there is no sudden un- 
bundling transition as in the models of Refs. 13, In 
those systems, it is the entropy stored in the transverse 
bending fluctuations that drives the transition towards 
the unbundled state. 

Let us now discuss the crosslink occupation as a func- 
tion of bundle deformation. For specifity, I take 8(s) = 
cN x sin(<7s) with q = n/2L, corresponding to an approx- 
imate solution for a clamped bundle that is loaded at 
its free end. The parameter c reflects the curvature of 



the bundle, 



and will serve in the following as a 



measure of the amplitude of the imposed deformation. 

As is evident from Fig. [2p,c the crosslink occupation 
(n) decreases upon increasing bundle deformation c. This 
agrees well with the above mentioned notion of bending- 
induced mismatch between the binding sites. An increas- 
ing mismatch implies increasing strain on the crosslink, 
which, in turn, increases their tendency to unbind. Un- 
like the case discussed in Fig.[2k, however, the decrease of 
(n) is gradual and smooth only below a certain value of 
the crosslink stiffness k x . Above this value, (n) displays 
a very sudden and discontinuous drop with bundle defor- 
mation. Interestingly, this effect cannot be explained on 
the basis of the one-crosslink toy- model of Eq. ((4]) , which 
only gives a smooth decay governed by a renormalized 
chemical potential 
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One concludes, that the discontinuous transition repre- 
sents a truly cooperative effect, where many crosslink 
unbinding events are needed to drive the bundle into a 
weakly crosslinked state. 

Analytical progress can be made in a mean-field ap- 
proximation. To this end we go back to the origi- 
nal Hamiltonian of Eq. ([2]) and substitute the actual 
crosslink state ni a by its bundle averaged value, n = 
J2ia n ia/(Nf — 1)N X . As compared to the fully occupied 
bundle (n = 1), this has the effect of renormalizing the 
crosslink stiffness k x — > nk x . 



The advantage of this approximation is that the it- 
variables can now be integrated out by going to the con- 
tinuum limit and transforming to Fourier space. The 
resulting effective free energy (per crosslink), which still 
depends on the bundle averaged occupation variable n 
has the rather simple structure [22[ 



F e g{n) = 



fj,n 



1 + 7i /n 

+ k B T [n log(n) + (1 - n) log(l - n)] , 



(7) 



Next to the chemical potential fi, there are two rele- 
vant parameters. The crosslink stiffness is encoded in 
n ~ (k s /k x )N f(N f + 1)/N X , while the parameter A 
is related to the amplitude of the bending deformation, 
A ~ k s b 2 c 2 Nf(N f + l). The term with the logarithm cor- 
responds to the standard entropy of mixing and counts 
the number of states that are compatible with a given av- 
erage crosslink occupation. Calculating the average (n) 
from Eq. ([7]) can easily be performed numerically, and 
is plotted as solid lines in Fig|2p,c. As can be seen the 
model and the simulation are in good agreement without 
fit parameters. 

At large bending amplitudes some deviations are ap- 
parent. By looking more closely into the distribution of 
crosslinks inside the bundle one finds that this state is ac- 
tually highly inhomogeneous, with many more crosslinks 
bound close to the clamped end than at the free end. 
In fact, the discontinuous reduction of bound crosslinks 
proceeds (with our MC-" dynamics" ) via an unzipping 
mechanism, where the crosslinks at the free end unbind 
first. An unbinding front then moves rapidly towards the 
clamped end of the bundle where it stabilizes some dis- 
tance away. A similar inhomogeneity can occur in the 
lateral direction, with crosslinks on different filaments 
unbinding at different bending amplitudes, leading to a 
sequence of discontinuous transitions, each correspond- 
ing to the unbinding of a single filament pair. I leave the 
detailed discussion of these complex unbinding pathways 
for a future publication. 

More insight into the physical mechanism that leads 
to the observed discontinuous transition can be obtained 
with the help of a saddle-point analysis. By minimizing 
Eq. ([7]) with respect to n, dF e s/dn = 0, one obtains the 
phase diagram depicted in Fig. |3] At T = (left figure) , 
i.e. without the entrop of mixing, the model has two 
minima at (n) =0,1 that coexist in a certain range of 
parameter values (limited by the dashed black lines). The 
condition F c g (0)=.F c ff (1) then gives a line of discontinu- 
ous transitions (solid red line) where the crosslink occu- 
pation jumps from (n) = 1 to (n) = 0. For large crosslink 
stiffness, the bending amplitude at the transition depends 
on the number of filaments as c x ~ (Nf(Nf + l)) -1 / 2 
(vertical dashed line). 

At finite temperatures (right figure) this line termi- 
nates at a critical point, beyond which the crosslink oc- 
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FIG. 3: Phase-diagram from a saddle-point analysis of Eq. ||?J| 
at T — (left) and for finite temperatures (right). In the 
latter case, the discontinuous transition (solid red line) ter- 
minates at a critical point, which shifts with the chemical po- 
tential along the (dotted) blue line. Metastable states exist 
within region delineated by (dashed) black lines. In the limit 
of large crosslink stiffness fc x — ► oo (thin vertical line) the dis- 
continuous transition asymptotes at Coo ~ l/y/Nf{Nf + 1). 



cupation decreases smoothly with increasing bending am- 
plitude. The critical point itself depends on the chemical 
potential [i as indicated in the figure by the (dotted) blue 
line. Thus, stiff crosslinks tend to unbind cooperatively 
in a discontinuous transition, while soft crosslinks unbind 
one after the other leading to a smooth decrease of the 
average crosslink occupation. 

This qualitative difference between bundles with soft 
or stiff crosslinks can be intuitively understood as fol- 
lows. When the crosslinks are soft, they cannot induce 
any stretching in the filaments. Accordingly the force 
in the crosslinks is small and set by the crosslink stiff- 
ness fc x - In this limit unbinding events are purely local 
with no consequences for the rest of the bundle. In con- 
trast, stiff crosslinks are able to drive the filaments into 
a stretched configuration. The force in the crosslink is 
now set by the filament stretching stiffness k s . An un- 
binding event will then affect the force balance in the 
filaments, with the potential of influencing the bundle 
state far away from the unbinding site. Note, that a sim- 
ilar distinction between soft and stiff crosslinks has been 



observed in Ref. [13], dealing with the phase behavior of 
helical filament bundles. 

The discontinuous transition and the associated hys- 
teresis observed in the simulations (see Fig|2|)b) is tan- 
tamount to a long time-scale associated with the es- 
cape over a free-energy barrier that separates highly 
crosslinked from weakly crosslinked states. Thus, the 
choice of crosslinker not only determines the mechanical 
properties of a bundle but also affects the relevant relax- 
ation times. Depending on the stiffness of the crosslink, 
this time-scale may range from the single crosslink bind- 
ing rate to the escape time over the free-energy barrier. 

Such a long time-scale may find its analogy in a re- 
cent experiment with F-actin bundles crosslinked by a- 



actinin [l8(. In the experiment a bundle was released 
from a bent configuration after keeping it there for a cer- 
tain waiting time. If the waiting time was long enough, 
the bundle did not relax back into the expected straight 
state, but remained with a considerable residual bend- 
ing deformation. It seems that under bundle deformation 
new binding sites become available that stabilize the bent 
shape by allowing the crosslinks to rebind in a "native" , 
unstrained state. Interestingly, the apparent time-scale 
(the waiting time) necessary to observe residual bundle 
bending was considerably longer than the time required 
for (un)binding of the individual cv-actinin linker. By 
taking /i = one can set up a simple correspondence be- 
tween unbinding events as defined in the present work, 
and rebinding events into the new sites in the experiment. 
With this mapping the long time-scale observed in the ex- 
periment can be identified with the escape time over the 
free energy barrier, which is estimated as r m 2 to * /N tot , 
where -/Vtnt = (Nf — 1)AT X is the total number of bind- 
ing sites [ISl, |20(. In this spirit, one would call a-actinin, 
the crosslinker used in the experiment, a stiff crosslink. 
It would be interesting to repeat the experiment with a 
softer crosslink. Future experiments should furthermore 
measure the force needed to bend the bundle, as this 
is expected to scale with crosslink occupation and dis- 
play a discontinuous drop, AF ~ N'j [201 ] . In particular, 
it would be interesting to discuss rate-dependent bundle 
bending. This would be a natural next step in the effort 
to understand the full dynamics of the coupled processes 
of bundle conformation and crosslink binding. Finally, 
it should be pointed out that crosslinked actin networks 
typically consist of composite mixtures of single filaments 
and small bundles Q . Therefore, the observed unbinding 
transition could directly be relevant for F-actin network 
rhcological properties. 
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